Quantum Tunneling Effect in Oscillating 
Friedmann Cosmology 

Mariusz P. Dabrowski 1 and Arne L. Larsen 2 
February 7, 2008 



Institute of Physics, University of Szczecin, Wielkopolska 15, 70-451 Szczecin 
Poland. 

2 Observatoire de Paris, DEMIRM. Laboratoire Associe au CNRS UA 336, Ob 
servatoire de Paris et Ecole Normale Superieure. 61, Avenue de l'Observatoire 
75014 Paris, France. 



Abstract 



It is shown that the tunneling effect in quantum cosmology is possible not 
only at the very beginning or the very end of the evolution, but also at 
the moment of maximum expansion of the universe. A positive curvature ex- 
panding Friedmann universe changes its state of evolution spontaneously and 
completely, without any changes in the matter content, avoiding recollapse, 
and falling into oscillations between the nonzero values of the scale factor. 
On the other hand, an oscillating nonsingular universe can tunnel sponta- 
neously to a recollapsing regime. The probability of such kind of tunneling 
is given explicitly. It is inversely related to the amount of nonrelativistic 
matter (dust), and grows from a certain fixed value to unity if the negative 
cosmological constant approaches zero. 



PACS number(s): 04.20.Jb; 04.60.+n; 98.80.Cq; 98.80.Hw 



1 Introduction 



The idea of unification of all forces in nature is quite old. It seems that GUT's 
or, in particular, the strong and electroweak theories have been accepted as 
standard theories of physics. However, the most interesting problem today 
is the unification of quantum mechanics with Einstein's theory of gravity. 
Some fundamentals of such a unification were first developed in the 60's 
by DeWitt, Wheeler and Misner U, @, |]. The idea that the universe could 
have been created through a quantum mechanical tunneling process from the 
vacuum, appeared in the 70's ||, |5], ||. More interest was concentrated on 
the topic in the 80 's, mainly because of the emergence of the idea of inflation 
[0] . Especially, Atkatz and Pagels considered a possibility for the universe 
to tunnel quantum mechanically from an initial static state, while Vilenkin 
U considered the tunneling from literally "nothing". Almost simultaneously, 
Hartle and Hawking |10| developed a formalism of quantum gravity which is 
based on the Feynman path integral approach. 

Recently, Dabrowski [[LlJ considered a class of Friedmann cosmologies 



which are oscillating, i.e., in which the universe oscillates between the nonzero 
values of the scale factor for infinitely long time. These solutions, in a way, 
generalize the static Einstein solution since they require some balance be- 



tween gravitational attraction and repulsion. In his original paper [12], Ein 



stein admitted nonrelativistic matter as a source of attraction and a positive 
cosmological constant as a source of repulsion. In our approach the prob- 
lem is more complex since we admit more sources corresponding to either 
attraction or repulsion. These are: nonrelativistic matter (attraction), the 
negative cosmological constant (attraction) and domain-wall-like matter (re- 
pulsion). In general, we can admit also radiation pressure (attraction), the 
positive cosmological constant (repulsion) and string-like-matter (repulsion), 
but they do not lead to any qualitative changes for the oscillating solutions 



II]. It should be stressed that the presence of do main- wall-like matter is es- 
sential for the existence of the nonsingular oscillating universes considered in 
this paper. Being quite a good candidate for the dark matter, stable domain 
walls are often considered to be cosmologically "disastrous" However, 
Hill et al [TJj] introduced the socalled "light" or "soft" (and thick) domain 
walls in a post-decoupling phase transition scheme, and found that the pro- 
duced domain walls were not necessarily in contradiction with the observed 
large scale structure of the universe (see also Ref . |L5], [l6j ) . We also point out 
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that strings and walls could have played a more important role in the early 
universe, where the tunneling processess under consideration here could be 
relevant, than they do today. Notice also that there is a way to avoid the do- 



main walls by replacing them by a kinematically equivalent scalar field [ 1 I 



which is postulated to be another source of the dark matter in the present era 
of the universe |17], For more discusion on the string-like and wall-like 
matter, see Ref. [11 

Usually, when dealing with quantum tunneling of the universe, all the 
quantum effects are considered to appear either at the very beginning or at 
the very end of the evolution. Recently, Kiefer and Zeh [|19| supported an 



earlier idea of Gold pOfl , that quantum effects for the universe as a whole 
may also be very important at the moment of its maximum expansion (i.e., 
when the universe changes its dynamics from expansion to collapse). Our 
work follows the spirit of these references and is most naturally related to the 
early paper by Atkats and Pagels ||, in the sense that we do not consider 
the quantum creation of the universe as a whole out of "nothing", but the 
quantum tunneling of the universe from one state to another. An important 
result of our analysis is that the evolution of the universe can change char- 
acter completely without any changes in the matter content. This should 
be contrasted with ordinary cosmological models where dramatic continuous 
(or discontinuous) changes in the scale factor usually follow from dramatic 
continuous (or discontinuous) changes in the matter content. Notice, in par- 



ticular, the paper by Hawking and Moss pi] , where a tunneling of the scalar 



field, equivalent to a discontinuous change in the matter content, leads to a 
spontaneous discontinous change of the Hubble parameter, i.e., of the evolu- 
tion of the universe. 

The plan of the paper is as follows. In Section 2 we comment on the clas- 



sically oscillating Friedmann cosmologies discussed previously in Ref. [[11 
Then we give a particular solution in terms of Jacobian Elliptic Functions, 
useful for the further purposes of this work. In Section 3 we use WKB ap- 
proximation in order to calculate the tunneling probability of the Friedmann 
universe from the expanding regime to the oscillating regime and from the 
oscillating regime to the recollapsing regime, for the particular solution given 
in Section 2. In Section 4 we discuss the results. 
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2 Classically Oscillating Friedmann Cosmolo- 
gies 

The existence of oscillating nonsingular solutions of the Friedmann equa- 



tion was first mentioned by Harrison [22]. He investigated qualitatively the 



universes filled with different kinds of both negative and positive pressure 
matter following the early discussion by Robertson |||] in the 30's. Recently, 
oscillating Friedmann cosmologies have been discussed qualitatively in the 



context of "exotic" matter by Kardashev |[24|| . A complete discussion of the 
oscillating solutions in terms of the Weierstrass Elliptic Functions has been 
given by Dabrowski fll|| . In this section we will consider in more detail a 



family of the simplest, but still nontrivial, solutions found in there. To fix 
our notation, and for later use, we first review a few general aspects of the 
construction. 

The nonvanishing components of the Ricci tensor for the generic four- 
dimensional Friedmann- Robertson- Walker line element, in comoving coordi- 
nates, 

dxdx 

ds 2 = -dt 2 + a 2 (t) (2.1) 
(1 + ^xx) 2 



a tt „ 2K + aa ttt + 2a 2 



are given by 

Rtt = -3—, Ru = — k ->->\2 ' "> ( 2 - 2 ) 
a (1 + ijxx) 1 

where a = a(t) is the scale factor, a t t = da/dt, a ttt = d 2 a/dt 2 and K is the 
curvature index. The corresponding scalar curvature is 

R= \(aa tt + a 2 t + K). (2.3) 

Including a nonvanishing cosmological constant A, the action takes the form 

S = S g + S m = -^— ( ^(R-2A)dQ + S m , (2.4) 

and the matter energy-momentum tensor Tn U is defined through 



5S, n 



IjV^g T^Sg^dQ. (2.5) 
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In Section 3, we shall return to the explicit expression for S m ; at this point 
we need only eq.(2.5), which together with the variation of eq.(2.4), leads 
to the Einstein equations in standard form. In the presence of radiation, 
nonrelativistic matter, string-like and wall-like matter, with corresponding 
densities C r , C m , C s and C w , respectively, the (tt)-component of the Einstein 



equations becomes [11, 25 




a 2 a 2 t = C r + aC m - a 2 (K - C s ) + a 3 C w + K\ (2.6) 

and the (ii)-component follows by differentiation with respect to t. Introduc- 
ing the conformal time, dr = dt/a(t), eq.(2.6) can be written as 

= «T 4 + \t 3 -{K- C s )T 2 + /3T + -, (2.7) 
where 

h~ c ^ = ho m , (3 = C W X- C 1 ' 2 , a = C r A c , A = A/A C , (2.8) 

and the so-called reduced temperature T is defined by 

A-i/2 

T=^ — , (2.9) 
a 

in agreement with Ref.]TT[. The Einstein equation is now of the elliptic 



form and can be explicitly solved in the general case. As shown in Ref.|Tl 
oscillating nonsingular solutions to eq.(2.7) can be obtained for a variety of 
parameter values. In the following we will consider the special, but nontrivial 

case 

a = 0, K=l, C s = 0, (3 = 3/8, (2.10) 

where the mathematics simplifies considerably; our results can however be 
easily extended to the more general case. We are thus considering a spatially 
closed universe with "ordinary" nonrelativistic matter and "exotic" wall-like 
matter. Equation (2.7) reduces to 

_ 2 T 3_ T2 3 A 
3 8 3' V ; 
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with the general solution 



Ar 1 / 2 



or 



6p(r - r ) + 1/2 



Here p is the Weierstrass Elliptic Function [EH] with invariants 



92 



1 

48' 



1 



93 



A 



1728 108 



and the discriminant reads as 



(2.12) 



(2.13) 



A^-27, 2 = -(A + -). 



In the special case A = —1/8, eq.(2.11) factorizes to [l 1 



(2.14) 



(2.15) 



and A = 0. In this case the solution (2.12) reduces to an expression in 
elementary functions PE| . Oscillating nonsingular solutions are obtained for 
A G ] — 1/8, 0[ , in which case A > 0. The three roots are given by 



ei 



1 + Z 2 



24 Z 



e 2 



iV3) (l + iV$)Z 



48 Z 



48 



e 3 



i\fl)Z 



48 Z 



where 



Z = Z(X) = [-1 - 16A + 4y/2^\{l + 8A) ] 1/3 . 
In the interesting region, A G ] — 1/8, 0[ , they fulfill 

ei > e 2 > e 3 . 



48 

(2.16) 
(2.17) 



(2.18) 



It is now straightforward to write the solution (2.12) in terms of Jacobian 
Elliptic Functions f2|, Hi 



or 



A c 1/2 sn 2 [yei - e 3 (r - r ) | (e 2 - e 3 )/(ei - e 3 )] 

6(ei - e 3 ) + (6e 3 + l/2)sn 2 [ v /e7 Tr e 3 "(r - r ) | (e 2 - e 3 )/ (e 1 - e 3 )] 

(2.19) 
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The integration constant tq must be carefully chosen to give real solutions 
for real r. Up to real translations, there are essentially two possibilities. For 
t = we find the solution 

/ v = A~ 1/2 sn 2 [yei - e 3 t | (e 2 - e 3 )/{e 1 - e 3 )] 

a{ ~ )[T) ~ 6( ei - e 3 ) + (6e 3 + l/2)s^[^^ r | (e 2 - e 3 )/( ei - e 3 )] ' 

(2.20) 

which oscillates between a = and a = Oj, where 

A-V2 

Taking instead r = iK'[(e 2 — e 3 )/(ei — e 3 )]/-^/ei — e 3 , where i^' is the 
complete elliptic integral of the first kind [26|, we get the solution 

A-i/2 

(+H J ~ (6e 3 + 1/2) + 6(e 2 - e 3 )sn 2 [ v ^7^ 3 " r | (e 2 - e 3 )/(e 1 - e 3 )] ' 

(2.22) 

which oscillates between a = a 2 and a = a 3 , where 

A -i/2 A -i/2 

" 2 ^ a ' s SETTS' (2 ' 23) 

The physical interpretation of these two solutions, a(_)(r), a( + )(r), follows 
by returning to the original equation of motion (2.6). For the parameter 
values (2.10), it takes the form 

( da\ 2 . i . i —2 3\/A^ a A c Aa 2 

UJ +K(o) = 0i K(0) = v^ +1 -^-^- (2 ' 24) 

Defining the potential in this way, we obtain that the dynamics takes place 
at the a-axis in a (a, l / (a))-diagram. For A G ] — 1/8, 0[ , the potential is 
shown in Fig.l, which explains the two solutions we found. In this picture, 
the universe is inflating (da/dt > and d 2 a/dt 2 > 0) when it is expanding 
and dV/da < 0. The solution a(_)(r) thus describes a deflating universe 
with scale factor expanding from a = to the maximal value a = a\, and 
thereafter re-collapsing to a = 0. The solution a( + )(r), on the other hand, 
describes a nonsingular universe, with scale factor oscillating between a = a 2 
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and a = a 3 . Notice also the behaviour in the two limits 



A — > : a% — > 0,2, as — > 00, 

A -> -1/8 : a 2 -> a 3 . (2.25) 

The comoving time of the solutions are obtained from 

t = f T a(T')dr', (2.26) 



which leads to 

, ( , K 1 ' 2 f 1 m 6e 3 + l/2 1 e a -e 3 l l 

6e 3 + 1/2 [ Vei - e 3 6(e 3 - d) e t - e 3 J 

(2.27) 

* ^ A ^ 1/2 rrr 6(e 3 -e 2 ) - /- — e 2 - e 3 

*(+)( r ) = 77 , ; /0 x ; =li « rTTo"' v e i ~ e 3 r > I 2 - 28 ) 

(6e 3 + l/2)y/ei - e 3 6e 3 + 1/2 e x - e 3 

where II is the elliptic integral of the third kind, and we are using the notation 
of Abramowitz and Stegun [2(J. For both solutions, the period in conformal 
time is given by 

P T = J 2 K[^^). (2.29) 
y /e 1 - e 3 ei - e 3 

It is however more relevant to consider the periods in comoving time t 

2A~V* l Kl e *- e h m 6e 3 + l/ 2 ■ e 2 -e 3 \ 

() (6e 3 + l/2) v ^7^ 3 " 1 [ e!-e 3 J [ 6(e 3 - e x ) 1 e x - e 3 J J ' 

(2.30) 

P (+) = , 1 fiZ*]. (2 . 31 ) 

(6e 3 + l/2) v ^T r ^ '6e 3 + 1/2 1 e t - e 3 J 

In the two extreme limits (2.25) we find 



A — > : P(_) — >■ oo, P( + ) — > oo, 
A -> -1/8 : P ( _) -► 47^(2 - v / 3)A f T 1/2 , P (+) -> Sttv^A,: 172 (2.32) 

It is interesting that the period Pu-\ is not continously going to zero for 
A — ► —1/8, as could have been expected from the potential in that particular 
limit (Fig.l), where the a(+) solution approaches the Einstein static Universe. 
This concludes our discussion of the classical behaviour of the oscillating 
solutions, eqs.(2.20),(2.22). 
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3 The Quantum Tunneling Probability 



As it is suggested from the potential given in Fig.l, although classically 
stable, the oscillating nonsingular solution a( + ) has the possibility to tunnel 
quantum mechanically through the barrier from a<i to a\ and then to collapse 
into a = 0. Similarly, the expanding solution a(_) which classically recollapses 
after hitting the barrier, has the possibility to tunnel quantum mechanically 
through the barrier from a\ to 02 and become oscillating and nonsingular. 
We believe that such processess could be relevant at least for the early evo- 
lution of the universe, i.e., shortly after the big bang. Notice that we are 
not discussing here the quantum creation of the universe as a whole out of 
"nothing" ||. Our aim is to consider a spontaneous change of character of 
evolution due to quantum processess after the big bang. In the standard 
cosmological models, the evolution of the universe changes because the mat- 
ter content changes, essentially because the temperature changes. What we 
want to illustrate with our model in this section is that there may be an 
additional effect: without any changes in the matter content, the evolution 
of the universe can change completely and spontaneously due to a quantum 
mechanical tunneling process. To be more specific, we will now calculate the 
probability for the universe to tunnel through the barrier between a = a\ and 
a = ci2, Fig.l. In the standard quantum mechanics calculation, the probabil- 
ity is given in the WKB-approximation by (see for instance the recent review 



by Atkatz, [28|, with applications to cosmology) 

p~e- B , (3.1) 
where B is twice the conjugate momentum integrated under the barrier 

(•0,2 

B = 2 \P a \da. (3.2) 

J a± 

Denoting by L, L g and L m , respectively, the Lagrangians corresponding to S, 
S g and S m , and assuming that the matter Lagrangian depends on the scale 
factor only via a, and not its derivatives (which is usually the case), we find 

_ dL dL g 
da >t da t t 

thus we do not need an explicit expression for the matter Lagrangian to 
calculate the tunneling probability (as a check, we calculate in Appendix 
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A the tunneling probability using the instanton method where S m is needed 
explicitly). The gravitational Lagrangian is obtained from eqs.(2.3),(2.4), for 
K = 1 

S g = ^L [ dt[-aa 2 t + a- Aa 3 /3], (3.4) 

4Cr J 

so that 

37T 

Pa = ~2G aa ' t ' 

For the parameter values (2.10), the Einstein equation (2.6) now takes the 
form 

?l = tM )2[ \ { ^ C a) ~ ^ + 3 8 { ^ C ^ + a)4] ' (3 ' 6) 

Then eq.(3.2) becomes 



/ \ rx% I 

B = J^q\]^~ J yx(x - xi)(x 2 - x)(x 3 - x) dx, (3.7) 

where Xi = y/A^. cij ; i — 1,2, 3. The integral is of elliptic type and the result 
can be expressed in terms of complete elliptic integrals |26|, |2{J, see Appendix 
B. For A — > 0, where a\ — > a 2 , we find as expected B = 0. There is no barrier 
in that limit so the "tunneling probability" is unity, but it actually takes 
infinite comoving time for the universe to reach the point a = a\ = a 2 (from 
below or above). In the other extreme limit, A — ► —1/8, where the solution 
0(+) is static, the integral (3.7) factorizes to give 



B(\ = -l/8) = —^-j A _ h (4-aA^a(a-A" 2 )da, (3.8) 
and we find 

From eq.(3.7) we conclude that the probability of tunneling (3.1), is inversely 
related to the overall amount of nonrelativistic matter (dust): for a small 
amount of nonrelativistic matter the tunneling probability approaches unity, 
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while it vanishes asymptotically for a large amount of it. This is quite reason- 
able since the quantum effects are thought to be connected with vacuum-like 
matter domination (cosmological constant, strings, walls). 

The plot of B as a function of A in the full range A G ] — 1/8, 0[, is shown 
in Fig. 2. For small (negative) A, the relationship is essentially linear while B 
increases steeply when A approaches —1/8. 



4 Conclusion 

In this paper we discussed the tunneling probability for a closed Friedmann 
expanding universe, to change its evolution completely and spontaneously 
and become an oscillating universe, without any changes in the matter con- 
tent. The tunneling happens merely when the closed universe reaches its 
maximum expansion size a\ (Fig.l). It falls into the oscillating regime, be- 
ginning its further evolution with a minimum value of the scale factor 02 
(ci2 > ai), and then oscillates between 02 and 03. Eventually, after some pe- 
riod of oscillations, when at the minimum 02 again, the universe can tunnel 
spontaneously without any changes in the matter content to the recollapsing 
regime. It begins collapsing exactly with the same value of the scale factor 
a\ corresponding to the maximum size of the classical Friedmann evolution. 
For the special solution of the Friedmann equation with fixed matter content 
of Section 2, we showed that the tunneling probability depends on the (here 
necessarily negative) cosmological constant and it is bigger for A closer to 
zero reaching unity for A = 0. The latter case corresponds to an unstable 
Einstein static solution together with the asymptotic solutions, see eq.(2.25) 
and Fig.l. 

Also, we found that the tunneling probability is inversely related to the 
amount of nonrelativistic matter (dust). It means that for a universe filled 
with more vacuum-type-matter the quantum effects are more probable and 
more important. This is a quite reasonable result, since the quantum effects 
are thought to be connected with vacuum-like matter domination (cosmolog- 
ical constant, strings, walls). 

Finally, one should mention the problem of arrow of time. Recently, 



Kiefer and Zeh [|19|, |30[ claimed that the arrow of time reverses just at the 
moment of maximum expansion a\ (Fig.l), in contradiction to the point of 
view of Hawking and collaborators (e.g. [[H], B3]). We do not discuss here the 
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question, whether there is a time asymmetry in our solutions, i.e., whether 
the expanding solution before the tunneling and the collapsing solution after 
the re-tunneling are described by the different arrow of time directions. It 
seems to be matter for further considerations. 
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A The Euclidean Action 



As was first shown by Coleman [33], the quantum mechanical tunneling prob- 



ability in the WKB-approximation, can also be obtained by an instanton 
calculation. The probability is 

p~e~ SE , (A.l) 

where Se is the Euclidean action of the classical solution in the classically 
forbidden region. In standard cases, the two expressions (3.1) and (A.l) lead 
to the same result. In the case of tunneling Friedmann universes, under con- 
sideration here, it is however instructive to perform the computation in both 
ways. The reason is that when using (3.1), we do not need an explicit expres- 
sion for S m (under the assumption that dS m /da jt = 0) and the computation 
proceeds as explained in Section 3 with the result (3.7). When using now 
(A.l) instead, we need an explicit expression for the matter action S m . This 
is usually not a problem when matter is represented by scalar (or higher spin) 
fields, but here the matter has been introduced only via energy and pressure 
densities in the energy-momentum tensor (2.5). An explicit expression for 
S m can still be obtained, but in a somewhat tricky way by introducing an 
additional function c(t) besides the scale factor a(t), see for instance Atkatz 
and Pagels pi. The action is written in the following way 

~ 37T f 1 2 n A , , 

9 = Tn I ~i aa ,t + ac ~ T7 a \ i 
4Gr J c 3 
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Sri 



3a/a: 



3tt r , o r 2 

•— — - / (it ca\ — ; o + 

AG J l 3^ c a 3 8a 



(A.2) 



S g reduces to eq.(3.4) for c = 1 and S m was chosen such that the Einstein 
equation (2.6), with the parameter values (2.10), appears when taking vari- 
ation with respect to c (and taking c = 1) in the total action S = S g + S m . 
Taking variation with respect to a, (and taking c = 1) we get simply the time 
derivative of the Einstein equation (2.6), with the parameter values (2.10). 
For c = 1, and defining t = it E , Se = iS, the Euclidean action is 



n 3-7T f r . da - o A o 2 

Se = 4G /period + ° " 3 ° " 



3v/A 



c 21 

a 



The Euclidean Einstein equation is 

-2a 



? 2/ \2 

(He 



A 



2 _ 3VA, a 3 _ ^ 4 
3a/A: 8 3 



(A.3) 



(A.4) 



Now it is straightforward to convert the t^-integral (A.3) into an a- integral, 
and via (A.l) we recover the result, eqs.(3.1),(3.7), for the tunneling proba- 
bility. 



B Exact Elliptic Integral, eq.(3.7) 

In this appendix, we give an explicit expression for the integral (3.7). Defining 



y 



x(x — Xx)(x2 — X)(X3 — x), 



and using the reduction formulas 17.1.4 in Ref . |[26| , we get 

3 r x * dx 1 81 r<n dx 

, L V 32^4 7~A^ ( 512 } 4 ~V 

1 2 27A 729 r** xdx 
+ A^ (A + "32~ + 8192 ) L ~V 



These are standard elliptic integrals 
1 r -3A 



ydx 



27A 729 



(B.l) 



(B.2) 



'XI 



X 3 Jx 2 (x 3 - x x ] 



16 



2*,(A' + — + 



, (A + )]AT(m) 

8192' x 3 v 512 ;J v ; 



12 



+ 



+ 



2(x 2 - x 3 



\ 3 ^Jx 2 (x 3 - Xi) 
2(x 2 - x 3 ) 



[A 2 + 



27A 729 



32 8192 J 



n 



,X\ — X 2 
X 3 — Xi 



^ + 81 ] n ^3(gl-g2) 

X 3 X2X 3 Jx 2 (x 3 - ^i) 512 ^2(^3-^1) 



:b.3) 



in the notation of Abramowitz [26|. The integral (3.7) is thus expressed 
in terms of complete elliptic integrals of first and third kind. The elliptic 
modulus is here given by 



m = k z 



x 3 (x 2 - Xi) 
x 2 (x 3 - Xx)' 



(B.4) 
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Figure Captions 



Fig.l. The potential V(a), denned in eq.(2.24), for a generic value of A in 
the range ] — 1/8, 0[ . The solution a(_) oscillates between a = and 
a — a\ (singular solution), while a(+) oscillates between a = a<i and a = a 3 
(nonsingular solution). 



Fig.2. The integral B, defined in eq.(3.7), plotted here as a function of 
A G ] — 1/8, 0[ . The tunneling probability is p ~ Exp (-£?). 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/gr-qc/9504025vl 



This figure "figl-2.png" is available in "png" format from: 
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